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The existence of optimal bang-bang controls for GMxB contracts* 


P. Azimzadeh^ P.A. Forsyth ^ 


Abstract 


A large collection of financial contracts offering gnaranteed minimnm benefits are often posed as 
control problems, in which at any point in the solution domain, a control is able to take any one of an 
uncountable number of values from the admissible set. Often, such contracts specify that the holder 
exert control at a finite number of deterministic times. The existence of an optimal bang-bang control, an 
optimal control taking on only a finite subset of values from the admissible set, is a common assumption 
in the literature. In this case, the numerical complexity of searching for an optimal control is considerably 
reduced. However, no rigorous treatment as to when an optimal bang-bang control exists is present in the 
literature. We provide the reader with a bang-bang principle from which the existence of such a control 
can be established for contracts satisfying some simple conditions. The bang-bang principle relies on the 
convexity and monotonicity of the solution and is developed using basic results in convex analysis and 
parabolic partial differential equations. We show that a guaranteed lifelong withdrawal benefit (GLWB) 
contract admits an optimal bang-bang control. In particular, we find that the holder of a GLWB can 
maximize a writer’s losses by only ever performing nonwithdrawal, withdrawal at exactly the contract 
rate, or full surrender. We demonstrate that the related guaranteed minimum withdrawal benefit contract 
is not convexity preserving, and hence does not satisfy the bang-bang principle other than in certain 
degenerate cases. 

Keywords: bang-bang controls, GMxB guarantees, convex optimization, optimal stochastic control 


1 Introduction 

1.1 Main results 

A large collection of financial contracts offering guaranteed minimum benefits (GMxBs) are often posed as 
control problems [3], in which the control is able to take any one of an uncountable number of values from the 
admissible set at each point in its domain. For example, a contract featuring regular withdrawals may allow 
holders to withdraw any portion of their account. In the following, we consider a control which maximizes 
losses for the writer of the contract, hereafter referred to as an optimal eontrol. 

A typical example is a guaranteed minimum withdrawal benefit (GMWB). If withdrawals are allowed at 
any time (i.e. “continuously”), then the pricing problem can be formulated as a singular control [18, 9, 14, 15] 
or an impulse control [7] problem. 

In practice, the contract usually specifies that the control can only be exercised at a finite number of 
deterministic exercise times to < ti < • ■ ■ < tw-i [3, 8]. The procedure for pricing such a contract using 
dynamic programming proceeds backwards from the expiry time t^ as follows: 
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1. Given the solution as t —>■ the solution as t —>■ is acquired by solving an initial value problem. 

2. The solution as t —> is then determined by applying an optimal control, which is found by considering 

a collection of optimization problems. 

If, for example, a finite difference method is used to solve the initial value problem from to an 
optimal control is determined by solving an optimization problem at each grid node, in order to advance the 
solution to t~. Continuing in this way, we determine the solution at the initial time. 

If there exists an optimal bang-bang control, an optimal control taking on only a finite subset of values 
from the admissible set, the numerical algorithm simplifies considerably. The existence of such a control 
is a common assumption in insurance applications [2, 19, 13], although no rigorous treatment is present in 
the literature. In this paper, we will also consider a weaker condition, a bang-bang principle. In this case, 
although an optimal control is not necessarily a finite subset of values from the admissible set, we will see 
that a control having this property can result in a large reduction in computational complexity. 

Our main result in this paper is the specification of sufficient conditions which can be used to guarantee 
the existence of an optimal bang-bang control. This result relies on the convexity and monotonicity of the 
solution and follows from a combination of basic results in convex analysis and parabolic partial differential 
equations (PDEs). We demonstrate our results on two common contracts in the GMxB family: 

• The guaranteed lifelong withdrawal benefit (GLWB) (a.k.a. guaranteed minimum lifelong withdrawal 
benefits (GMLWB)) admits an optimal bang-bang control. In particular, we prove that a holder can 
maximize the writer’s losses by only ever performing 

— nonwithdrawal, 

— withdrawal at the contract rate (i.e. never subject to a penalty), or 

— a full surrender (i.e. maximal withdrawal; may be subject to a penalty). 

• On the other hand, the guaranteed minimum withdrawal benefit (GMWB) is not necessarily convexity 
preserving, and does not satisfy a bang-bang principle other than in certain degenerate cases. 

In the event that it is not possible to determine an optimal control analytically, numerical methods are 
required. Standard techniques in optimization are not always applicable, since these methods cannot guarantee 
convergence to a global extremum. In particular, without a priori knowledge about the objective functions 
appearing in the family of optimization problems corresponding to optimal holder behavior at the exercise 
times, a numerical method needs to resort to a linear search over a discretization of the admissible set. 
Gonvergence to a desired tolerance is achieved by refining this partition [23]. Only with this approach can 
we be assured of a convergent algorithm. However, if an optimal bang-bang control exists, discretizing 
the control set becomes unnecessary. Theoretically, this simplifies convergence analysis. More importantly, 
in practice, this reduces the amount of work per local optimization problem, often the bottleneck of any 
numerical method. 

1.2 Insurance applications 

The GLWB is a response to a general reduction in the availability of defined benefit pension plans [6], allowing 
the buyer to replicate the security of such a plan via a substitute. The GLWB is bootstrapped via a lump sum 
payment wq to an insurer, which is invested in risky assets. We term this the investment account. Associated 
with the GLWB contract is the guaranteed withdrawal benefit account, referred to as the withdrawal benefit 
for brevity. This account is also initially set to wq. At a finite set of deterministic withdrawal times, the 
holder is entitled to withdraw a predetermined fraction of the withdrawal benefit (or any lesser amount), 
even if the investment account diminishes to zero. This predetermined fraction is referred to as the contract 
withdrawal rate. If holders wish to withdraw in excess of the contract withdrawal rate, they can do so upon 
the payment of a penalty. Typical GLWB contracts include penalty rates that are decreasing functions of 
time. 
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These contracts are often bundled with ratchets (a.k.a. step-ups), a contract feature that periodically 
increases the withdrawal benefit to the investment account, provided that the latter has grown larger than 
the former. Moreover, bonus (a.k.a. roll-up) provisions are also often present, in which the withdrawal benefit 
is increased if the holder does not withdraw at a given withdrawal time. Upon death, the holder’s estate 
receives the entirety of the investment account. We show that a holder can maximize the writer’s costs by 
only ever performing nonwithdrawal, withdrawal at exactly the contract rate, or surrendering the entirety of 
their account. Such a holder will never withdraw a nonzero amount strictly below the contract rate or perform 
a partial surrender. However, this result requires a special form for the penalty and lapsation functions, which 
is not universal in all contracts. Pricing GLWB contracts has previously been considered in [21, 13, 10, 1]. 

Much like the GLWB contract, a GMWB is composed of an investment account and withdrawal benefit 
initially set to wq, in which wq is a lump sum payment to an insurer. At a finite set of withdrawal times, the 
holder is entitled to withdraw up to a predetermined amount. Note that this amount is not a fraction of the 
withdrawal benefit, as in the GLWB, but rather a constant amount irrespective of the withdrawal benefit’s 
size. Furthermore, unlike the GLWB, the action of withdrawing decreases both the investment account and 
withdrawal benefit on a dollar-for-dollar basis. 

The GMWB promises to return at least the entire original investment, regardless of the performance of 
the underlying risky investment. The holder may withdraw more than the predetermined amount subject to 
a penalty. Upon death, the contract is simply transferred to the holder’s estate, and hence mortality risk 
need not be considered. Pricing GMWB contracts has been previously considered in [18, 9, 8, 14, 15[. 

1.3 Overview 

In §2, we introduce the GLWB and GMWB contracts. In §3, we generalize this to model a contract that can 
be controlled at finitely many times, a typical case in insurance practice (e.g. yearly or quarterly exercise). 
In §4, we develop sufficient conditions for the existence of an optimal bang-bang control and show that the 
GLWB satisfies these conditions. §5 discusses a numerical method for finding the cost of funding GLWB and 
GMWB contracts, demonstrating the bang-bang principle for the former and providing an example of where 
it fails for the latter. 


2 Guaranteed minimum benefits (GMxBs) 


We introduce mathematical models for the GLWB and GMWB contracts in this section. Since most GMxB 
contracts offer withdrawals on anniversary dates, to simplify notation, we restrict our attention to annual 
withdrawals occurring at 

^ = {0,l,...,N-l}. 

0 and N are referred to as the initial and expiry times, respectively (no withdrawal occurs at N). 

In order to ensure that the writer can, at least in theory, hedge a short position in a GMxB with no 
risk, we assume that the holder will employ a loss-maximizing strategy. That is, the holder will act so as to 
maximize the cost of funding the GMxB. This represents the worst-case hedging cost for the writer. This 
worst-case cost is a function of the holder’s investment account and withdrawal benefit. As such, we write 
X = {xi,X 2 ), where Xi is the value of the investment account and X 2 is the value of the withdrawal benefit. 
Both of these quantities are nonnegative. 

Let a denote the hedging fee, the rate continuously deducted from the investment account Xi (while xi 
is used to denote a particular value of the investment account, the capital symbol Xi is reserved for the 
corresponding stochastic process) to provide the premium for the contract. We assume that between exercise 
times, the investment account of the GMxBs follows geometric Brownian motion (GBM) as per 


dXi 


{pL — ot)dt + adZ 
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tracking the index Xi satisfying 


dXi 


^dt+ adZ 


where Z is a, Wiener process under the real-world measure. We assume that it is not possible to short the 
investment account Xi for fiduciary reasons [8], so that the obvious arbitrage opportunity is prohibited. 

The worst-case cost of a GMxB is posed as the solution to an initial value problem (IVP) specified by 
three conditions: 


1. the worst-case cost of funding the contract at the expiry time (posed as a Cauchy boundary condition; 
see, for example, (2.1) and (2.11)); 

2. the evolution of the worst-case cost across withdrawals (posed as a supremum over the holder’s actions, 
corresponding to the holder acting so as to maximize the writer’s losses; see, for example, (2.2) and 
( 2 . 12 )); 

3. the evolution of the worst-case cost between withdrawals (posed as a conditional expectation; see, for 
example, (2.3) and (2.13)). 

We begin by introducing the IVP for the GMWB before moving to the GLWB for ease of exposition. To 
distinguish the two contracts, we use the superscripts L and M to denote quantities that pertain to the 
GLWB and GMWB, respectively. In the following, we denote by E the expectation and by Z a Wiener 
process under the risk-neutral measure, that which renders the discounted index Xi into a martingale. For a 
function g whose domain is a subset of K, we use the notations g (t~) = limg^t g (s) and g (t"*") = lirng^t g (s) 
to denote the one-sided limits at t. 


2.1 Guaranteed minimum withdrawal benefit (GMWB) 

Since the GMWB is transferred to the holder’s estate upon death, mortality risk is not considered. The 
worst-case cost of funding a GMWB at time N (the expiry) is [9] 

(x) = max (xi, (1 — kn) X 2 ), 


corresponding to the greater of the entirety of the investment account or a full surrender at the penalty rate 
at the iVth anniversary, hm S [Oj !]• The worst-case cost of funding a GMWB at previous times is derived 
by a hedging argument in which the writer takes a position in the index Xi [8]. Equivalently, it is given by 
finding V (within the relevant space of functions; see Appendix A) such that (s.t.) 


V (x,N) = (x) 

V (x, n-) = sup [V (f^„ (A), n+) -k (A)] 

Ag[0,1] 


V (x, t) = E 


e- /:" (Vi ((n + 1)-) , (n + 1)") 

I Xi (n+) = xi 


on [ 0 ,oo)^ 

(2.1) 

on [0, 00)^ X .X 

(2.2) 

on [0, 00 )^ X {n,n-\- 1) Vn 

(2.3) 


where between exercise times 

—^ = (r — a) dt -h adZ. (2.4) 

Xi 

r is the risk-free rate, : [0,1] —M represents the cash flow from the writer to the holder, and : [0,1] —>■ 
[0,oo) represents the state of the contract postwithdrawal. The construction of and is outlined below. 
The holder is able to withdraw a fraction A S [0,1] of the withdrawal benefit at each exercise time. 

Intuitively, V {x,n~) and V (x,n+) can be thought of as the value of the contract “immediately before” 
and “immediately after” the exercise time n. 
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Let G ^ 0 denote the predetermined eontraet withdrawal amount associated with the GMWB so that 
G A X2 (a A 6 = min (a, b), aM b = max (a, b)) is the maximum the holder can withdraw without incurring a 
penalty (both A and V are understood to have lower operator precedence than the arithmetic operations). 
Consider the point {xi,X 2 ,n) with n & 3A. 

• The maximum a holder can withdraw without incurring a penalty is G A a; 2 . If the holder withdraws 
the amount Aa ;2 with Ax 2 G [0, G A X2], 

V (x,n“) = V{xi — \x 2 V 0, ^2 ~ ^X 2 , n'^) + \x 2 ■ (2-5) 

fM /M 


• Let € [O 5 1] denote the penalty rate at the nth anniversary. If the holder withdraws the amount Xx 2 
with Aa ;2 G (G A 0 : 2 , 2 : 2 ], 

V (x,n~) = V{xi — Xx 2 V 0, X 2 — Xx 2 , n'^) + Xx 2 — Kn (Aa ;2 — G). (2-6) 

fM JM 

Here, Aa ;2 G (G A X 2 , X 2 ) corresponds to a partial surrender and Aa :2 = X 2 (i.e. A = 1) corresponds to a 
full surrender. 


We can summarize (2.5) and (2.6) by taking 

Aa;2 


/“„(A)^ 


and 


if Xx2 G [0, G A X2] 
i G + (1 — K„) (Ax 2 — G) if Aa;2 G (G A X2, X2] 

(A) = {xi - XX2 V 0, (1 - A) 2 : 2 ). 


(2.7) 


It can be shown from (2.3) that the cost to fund the GMWB (between exercise times) satisfies^ [ 8 ] 


where 


dtV + CV = 0 on (0, 00 )^ X (n, n + 1) Vn 


C = + (r-a) Xid^, - r. 


( 2 . 8 ) 

(2.9) 


2.2 Guaranteed lifelong withdrawal benefit (GLWB) 

Let At (t) be the mortality rate at time t (i.e. A4 (t) dt is the fraction of the original holders who pass 
away in the interval [ti,t 2 ]), so that the survival probability at time t is 

Tl{t) = 1 — f M{s) ds. 

Jo 

We assume X4 is continuous and nonnegative, along with 77. (t) ^ 0 for all times t. We assume that mortality 
risk is diversifiable. Furthermore, we assume the existence of a time t* > 0 s.t. 77 (t*) = 0. That is, survival 
beyond t* is impossible (i.e. no holder lives forever). N is chosen s.t. TV ^ t* to ensure that all holders have 
passed away at the expiry of the contract. As is often the case in practice, we assume ratchets are prescribed 
to occur on a subset of the anniversary dates (e.g. triennially). 

As usual, we assume that the holder of a GLWB will employ a loss-maximizing strategy. Since TV was 
picked sufficiently large, the insurer has no obligations at the TVth anniversary and the worst-case cost of 
funding a GLWB at time TV is 

y(x) = 0. (2.10) 

^We discuss what it means for a function to satisfy this PDE in Appendix A. 
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As with the GMWB, the worst-case cost of funding a GLWB is derived by a hedging argument in which the 
writer takes a position in the index Xi [10]. Equivalently, it is given by finding V (within the relevant space 
of functions; see Appendix A) s.t. 


V (x,N) = (x) 

V (x, n-) = sup [E (f^„ (A), n+) -P (A)] 

AG[0,2] 

E(x,t) ((n+1)-) ,X 2 ,(n + l)-) 

/ n+1 ps 

e Jt (s) Xi (s) ds I Xi {n~^) = Xi 


on [0,oo)^ (2.11) 

on [0,oo)^ X (2.12) 


on [0, oo)^ X (n, n-I-1) Vn (2.13) 


where between exercise times, Xi is specified by (2.4). 

: [0,2] —)■ K represents the (mortality-adjusted [10]) cash flow from the writer to the holder and 
f^: [0,2] —)■ [0,oo) represents the state of the contract postwithdrawal. In particular, A = 0 corresponds to 
nonwithdrawal, A G (0,1] corresponds to withdrawal at or below the contract rate, and A G (1, 2] corresponds 
to a partial or full surrender. 


Remark 2.1. We remark that the admissible set of actions [0,2] is undesirably large (i.e. a continuum). We 
will apply the results established in §4 to show that an optimal strategy taking on values only from {0,1,2} 
exists. In other words, an eguivalent problem can be constructed by substituting the set {0,1,2} for the original 
[0,2] in the optimization problem (2.12). The resulting problem has smaller computational complexity than 
the original one (i.e. successive refinements of [ 0 , 2 ] need not be considered to attain convergence). 

The construction of f^ and is guided by the specification of the contract: 


• Let /3 denote the bonus rate, if the holder does not withdraw, the withdrawal account is amplified by 

1 + /3. 

• Let 5 denote the contract withdrawal rate; that is, 6 x 2 is the maximum a holder can withdraw without 
incurring a penalty. 

• Let Kn € [0, 1] denote the penalty rate at the nth anniversary, incurred if the holder withdraws above 
the contract withdrawal rate. 


• Let 


Then, 


and 


In = 


1 if a ratchet is prescribed to occur on the nth anniversary 
0 otherwise 


0 

ftn {X)=n{n)- {XSX2 


if A = 0 
if A e (0,1] 


6 x 2 -I- (A — 1) (1 — Kn) (a^i — 6 x 2 VO) if A G (1,2] 


{xi, X2 (1 -b ,d) V InXl) 


if A = 0 


W,n \ ( 2:1 — XSx 2 V 0, 0:2 V I„ [xi — X 6 X 2 ]) if A G (0,1] . 


(2-A)fx,„(l) if A G (1,2] 

It can be shown that the cost to fund the GLWB (between exercise times) satisfies [10] 


dtV -b CV + Mxi = 0 on (0, 00 )^ x (n, n -b 1) Vn 


(2.14) 


(2.15) 


(2.16) 


where £ is defined in (2.9). 
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3 General formulation 


We generalize now the above IVPs. Let 3^ = {to, • • ■ along with the order 0 = to < ■ ■ ■ < tN = T, m 

which T is referred to as the expiry time. Let be a convex subset of K’". The set of all actions a holder 
can perform at an exercise time is denoted by A„ C K™ , assumed to be nonempty, and referred to as an 
admissible set. For brevity, let 

v^.n (A) = P (fx,„ (A) , t+) + /x,„ (A) (3.1) 

where : A„ —>■ K and fx,n: A„ —)■ 12. We write v^.n (A) to stress that for each fixed (x, n), we consider an 
optimization problem in the variable A. The general problem is to find V satisfying the conditions 

P (x, r) = (/? (x) on 12 (3.2) 

V (x,t~) = supux,„ (A„) on 12 X (3.3) 

along with a condition specifying the evolution of V from 2+ to t~_i_i (see, for example, (2.3) and (2.13)). 

Remark 3.1. Convexity preservation, a property that helps establish the bang-bang principle, depends on 
each admissible set Kn being independent of the state of the contract, x. This is discussed in Remark 4.18. 

4 Control reduction 

Definition 4.1 (Optimal bang-bang control). V, a solution to the general IVP introduced in §5, is said to 
admit an optimal bang-bang control at time tn G 3^ whenever 

V (x,2“) = maxux,„ (^A„^ on 12, 
where An denotes a finite set independent of x. 

The above condition is inherently simpler than (3.3), in which there are no guarantees on the cardinality 
of A„. 

§4.2 develops Corollary 4.13, establishing sufficient conditions for the existence of an optimal bang-bang 
control. This result requires that the relevant solution V be convex and monotone (CM). Given a CM initial 
condition (3.2), we seek to ensure that V preserves the CM property at all previous times. §4.3 develops 
conditions on the functions / and f to ensure that the supremum (3.3) preserves the CM property. Similarly, 
§4.4 develops conditions on the dynamics of V (and hence the underlying stochastic process(es)) to ensure 
that the CM property is preserved between exercise times. 

For the remainder of this work, we use the shorthand Vn (x) = V (x,2+) and Vf (x) = V {x,t~). 

4.1 Preliminaries 

In an effort to remain self-contained, we provide the reader with several elementary (but useful) definitions. 
In practice, we consider only vector spaces over M and hence restrict our definitions to this case. 

Definition 4.2 (convex set). Let W be a vector space over M. X C W is a convex set if for all x,x' G X 
and 9 € (0,1), Ox -\- {1 — 9) x' € X. 

Definition 4.3 (convex function). Let X be a convex set and Y be a vector space over K equipped with a 
partial order ^y. h: X ^ Y is a convex function if for all x,x' G X and 9 G (0,1), 

h {Ox -I- (1 — 0) x') ^y 9h (x) -\- {1 — 9)h {x'). 

Definition 4.4 (extreme point). An extreme point of a convex set X is a point x G X which cannot be 
written x = Ox' -I- (1 — 0) x" for any 9 G (0,1) and x', x" £ X with x' x". 
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Definition 4.5 (convex polytope). Let Y he a topological vector space over K. P <ZY is a convex polytope if 
it is a compact convex set with finitely many extreme points. The extreme points of a convex polytope are 
referred to as its vertices. 

Definition 4.6 (monotone function). Let X and Y be sets equipped with partial orders and 
respectively, h: X ^Y is monotone if for all x,x' G X, x x' implies h{x) h{x'). 

Lemma 4.7. Let A be a convex set, and let B and C he vector spaces over M equipped with partial orders 
and respectively. If hi: A ^ B and / 12 : B ^ C are convex functions with /12 monotone, then h 2 o hi 
is a convex function. 

Remark 4.8. For the remainder of this work, we equip ffi™ with the order ^ defined as follows: */x, y G IR™, 
X ^ y whenever Xi ^ yi for all i. 

4.2 Bang-bang principle 

Consider a particular exercise time tn. Suppose the following: 

(Al) X V+ (x) is CM. 

(A2) For each fixed x G Vx,n (A„) is bounded above. 

Throughout this section, we consider a particular point y G f2 in order to establish our result pointwise. For 
the results below, we require the following propositions: 

(Bl) There exists a collection Vn (y) C 2^" s.t. P = A„ and each P £ Vn (y) is compact convex. 

(B2) For each P GVn (y), the restrictions A i-G fy.n\p (A) and A fy,n|p (A) are convex. 

(B3) Vn (y) is a finite collection of convex polytopes. 

Remark 4.9. (Bl) simply states that we can “cut up” the admissible set An into (possibly overlapping) 
compact convex sets. (B2) states that the restrictions of /y^„ and {y,n on each of these sets are convex 
functions of X. 

Lemma 4.10. Suppose (Al), (Bl), and (B2). For each P G P„ (y), the restriction A 1 —>■ Wy,n|p (A) is convex. 
Proof. The proof is by (3.1), (Al), (B2), and Lemma 4.7. □ 

Lemma 4.11. Suppose (Al), (A2), (Bl), and (B2). Let P G Vn (y)- Then, 

sup Vy^n (P) = sup Vy^n {E (P) ) 
where E (P) denotes the set of extreme points of P. 

Proof. Let w = Vy^n\p- Note that w (P) = Uy^„ (P), and hence no generality is lost in considering w. Lemma 
4.10 establishes the convexity of w. Naturally, supu> (P) exists (and hence sup w {E (P)) exists too) due to 
(A2). Finally, it is well known from elementary convex analysis that the supremum of a convex function on a 
compact convex set P lies on the extreme points of P, E (P). See [22, Chap. 32]. □ 

Theorem 4.12 (bang-bang principle). Suppose (Al), (A2), (Bl), and (B2). Then, 


supuy,„ (A„) = supuy,„ 


U 

PeP„{y) 


where E (P) denotes the set of extreme points of P. 


Proof. By (Bl), we have that A„ = IJpgp^(y) P. We can, w.l.o.g., assume that all members of Vn (y) are 
nonempty (otherwise, remove all empty sets), supay^n (A„) exists due to (A2). Since for each P £ Vn (y), 
supuy^n {P) = supuy^n {E (P)) (Lemma 4.11), two applications of Lemma B.l allow us to “commute” the 
supremum with the union to get 


supuy^n (A„) = supuy^n I |J P 

KPev„{y) 


= supay,„( IJ P(P)). 
KP^^'Pr.iy) 


□ 


Theorem 4.12 reduces the region over which to search for an optimal control. When Vn (y) is a finite 
collection of convex polytopes, the situation is even nicer, as Upgp (y) ^ ^ finite set (a finite union of 

finite sets). If, in addition, Vn is chosen independent of y, we arrive at an optimal bang-bang control: 

Corollary 4.13 (optimal bang-bang control). Suppose (Al) and (A2). Furthermore, suppose (Bl), (B2), 
and (B3) for all y G ft. Finally, suppose that there exists Vn s.t. Vn = Vn (y) for all y G ft. Then, the 
general IVP introduced in §5 admits an optimal bang-bang control at time tn (Definition A.l) with 

V (x,t“) = supux.n (A„) = maxUx.ri (a„^ on D. 

and 

A„= U E{P). 

PG'P„ 

Example 4.14. Let yG [0,oo)^. We now find V(^ (y) s.t. (Bl), (B2), and (B3) are satisfied for the GLWB. 
Take Pi = [0,1], P 2 = [1,2], and V(^ (y) = {Pi, P 2 }) satisfying (B3). Note that P = [0, 2], satisfying 

(Bl). It is trivial to show that the functions fyn\pj ond fy „|pj defined in (2.14) and (2.15) are convex as 
functions of X (the maximum of convex functions is a convex function), thereby satisfying (B2). Since y was 
arbitrary and V(i was chosen independent ofy, we conclude (whenever (Al) and (A2) hold), by Corollary 
4.13, that the supremum ofVy „ occurs at 

A(: = E (Pi) U E (P2) = P ([0,1]) U E ([1,2]) = (0,1} U (1, 2} = (0,1, 2} 

(corresponding to nonwithdrawal, withdrawal at exactly the contract rate, and a full surrender). 

Remark 4.15. When all the conditions required for Corollary 4.13 hold, with the exception that P„ (y) 
depends on y, then an optimal control is not necessarily bang-bang, but does satisfy the bang-bang principle. 
Theorem 4.12. In many cases, this still results in considerable computational simplification (see Remark 5.3). 


4.3 Preservation of convexity and monotonicity across exercise times 

Since the convexity and monotonicity of V are desirable properties upon which the bang-bang principle 
depends (i.e. (Al)), we would like to ensure that they are preserved “across” exercise times (i.e. from to 

if)- 

Consider the nth exercise time, tn. Suppose the following: 

(Cl) For each fixed A G A„, x i-> fx,„ (A) and x /x,„ (A) are convex.^. 

^Note that this is not the same as (B2) Here, we mean that for each fixed A S A„ and for all x,x' S O and 0 S (0,1), 

/ex+{l-e)x',n (-^l ^ (A) -I- (1 — 0) /x',n W 

and 

fex+(l-e)x',n W ^ 9tx,n (A) -I- (1 — 0) fx'_„ (A). (4.1) 

The order Si used in (4.1) is that on SI C R*", inherited from the order on R*" established in Remark 4.8. 
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(C2) For each x,x' e s.t. x ^ x', there exist sequences {Xk} s.t. Ux,n {Xk) —>■ (x), and for 

all k, /x,„ (Afc) < /x',„ (A'j.) and fx,„ (Afe) fx'.n (A^). 

Remark 4.16. (C2) simplifies greatly if for all x, Wx,n (A-n) contains its supremumf' Denote this supremum 
Vyi.n (Ax): wherc Ax G A„ is an optimal action at x. In this case, the following simpler assumption yields 
(C2); for each x,x' G s.t. x ^ x', there exists X' G A„ s.t. fy^^n (Ax) ^ fx',n (A') and fx,™ (Ax) ^ fx',n (A') 
(take Xk = Ax and A(. = A' for all k to arrive at (C2 ) ). 

This simpler condition states that for each pair of positions x ^ x', there is an action X' s.t. the position 
and cash flow after the event at x' under action A' are greater than (or equal to) the position and cash flow 
after the event at x under an optimal action Ax- Intuitively, this guarantees us that the position x' is more 
desirable than x (from the holder’s perspective). This is not a particularly restrictive assumption, and it 
should hold true for any model of a contract in which a larger position is more desirable than a smaller one. 

Lemma 4.17. Suppose (Al), (A2), and (Cl). Then, x i-G V~ (x) is convex. 

Proof. Fix X, x' G n and 9 G (0,1), and let z = 0x + (1 — 0) x'. Then, by (Al) and (Cl), 

V~ (z) = SUpUz,n (A„) 

= sup [V+ (fz.„ (A)) + /z,„ (A)] 

aga„ 

< sup [V+ (0fx.„ (A) + (1 - 0) fx'.n (A)) + 0/x.„ (A) + (1 - 0) f^>,n (A)] 

AgA„ 

^ 0 sup [V+ (fx.n (A)) + /x.n (A)] + (1 - 0) SUp [V+ (fx',„ (A)) + /x'.n (A)] 

AGA„ AgA„ 

= 0 sup Ux.n (A„ ) + (1 - 0) sup Ux' ,„ (A„ ) 

= ev- (x) + (1 - 0 ) y- (x'). 


□ 

Remark 4.18. Note that the proof of Lemma 4.17 involves using V~ (y) = supuy.„ (A„) for y = x, x'. If 
An is instead a function of the contract state (i.e. An = An (x )), then the above proof methodology does not 
work since it is not necessarily true that Vfl (y) = supuy.n (A„ (z)) for y = x,x'. 

Lemma 4.19. Suppose (Al), (A2), and (C2). Then, x i—> Vn (x) is monotone. 

Proof. Let x,x' G Al s.t. x ^ x'. By (Al) (specifically, since is monotone) and (C2), for each k, 

Vxi.n (Afc) = Vn (fx ,n (Afc)) + /x ,n (Afc) 

^ v;+ (fx'.n (A',)) + /x'.n (A'fc) 

= Vx'.n (Afc) 


Then, 


as desired. 


Vn (x) = lim Vx,n (Afe) < limsupux'.n (A'j.) ^ supux'.n (A„) = Vn (x') , 

k^oo fe->.oo 


□ 


Example 4.20. We now show that the GLWB satisfies (Cl) and (C2) given (Al) and (A2). It is trivial 
to show that the functions f^n (''^) i^n ('^) defined in (2.14) and (2.15) are convex in x (the maximum 
of convex functions is a convex function), thereby satisfying (Cl). (C2) is slightly more tedious to verify. 
Let x,x' G s.t. X ^ x'. By (Al), (A2) and the argument in Example 4.14, we can, w.l.o.g., assume 
Xx G {0; 1) 2}, where Ax denotes an optimal action at x. Hence, we need only consider three cases: 

®It is worthwhile to note that in practice, this is often the case; for fixed n, consider A„ compact and A i—>■ Vx.n (A) continuous 
for all X. 
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1 . Suppose Ax = 0. Take X = 0 to get (0) = „ (A') and f^„ (0) f^, „ (A')- 

2. Suppose Ax = 1. W.l.o.g., we can assume X2 ^ X2 > 0. Take X = x^jx'^ to get (1) = f^, ^ (A') and 

S. Suppose Ax = 2. If xi X, 5x2, then = /x.n(l) fx,n (2) = (0,0) ^ (1), and we 

can w.l.o.g. assume X2 ^ X2 > 0 and once again take X = X2lx'2 to get „ (2) = /x'„ (^0 
fx n (2) = (0, 0) ^ f^/ „ (A'). Therefore, we can safely assume that xi > 8x2 so that 

fx,n (2) = {n) [(1 - k)xi + k 5 x 2] < Tl (n) xi. (4.2) 

(a) Suppose x'l ^ 6x2. Take X = 1 to get (2) = (0,0) ^ „ (1) and 

ftn (2) X,'R{n)xiif'R (n) 8x2 = ft',n (1) 

by (4.2). 

(b) Suppose x'l > 5 x' 2. Take X = 2 to get ^ (2) = (0,0) = „ (2) and 


fLu (2) {n) [(1 - Kn) x'l + k5x'2\ = /x,,„ (2). 


4.4 Preservation of convexity and monotonicity between exercise times 

As previously mentioned, to apply Theorem 4.12, we need to check the validity of (Al) (i.e. that the solution 
is CM at t+). In light of this, we would like to identify scenarios in which Vff is CM provided that VX+i is 
CM (i.e. convexity and monotonicity are preserved between exercise times). 

Example 4.21. If we assume that both GLWB and GMWB are written on an asset that follows GBM, then 
Appendix A establishes the convexity and monotonicity (under sufficient regularity) ofVff given the convexity 
and monotonicity of V~_^i. The general argument is applicable to contracts written on assets whose returns 
follow multidimensional drift-diffusions with parameters independent of the level of the asset (a local volatility 
model, for example, is not included in this class). Convexity and monotonicity preservation are retrieved 
directly from a property of the corresponding Green’s function. 

Although the methodology in Appendix A relates convexity and monotonicity to a general property of the 
Green’s function (including the class of contracts driven by GBM), in the interest of intuition, we provide the 
reader with an alternate proof below using the linearity of the expectation operator along with the linearity of 
the stochastic process w.r.t. its initial value. Consider, in particular, the GLWB. Equation (2.13) stipulates 


V+ (x) = E 


e 


+ e~E (s) Xi (s) ds I Xi (n+) 

J n 


Xi 


on [0, c»)^ X . 


Linearity allows us to consider the two terms appearing in the .sum inside the conditional expectation separately. 
If each is convex in x, so too is the entire expression. If Xi (n^) = xi, 

Xi{s)=XiYis) 


between n and n + 1, where 


Y (s) = exp 


r(r) - a(T) - -cr^ (r) 


dr + / (T (r) dZ (r) ) , 


from which it is evident that Xi (s) is convex in xi since Y depends only on time (note that the parameters 
appearing in Y are independent of the level of the asset, precluding a local volatility model). It remains to 
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show that the first term is also convex. Fix y,y' G [0, oo)^, 9 G (0,1), and let x = 0y + (1 — 0) y'. Then, by 
assuming that (x) is convex in x, 

(xiF ((n+ 1)") ,^ 2 ) = Ki+i [{ 6 yi + {l-e)y[)Y ((n+ 1)") , 9 y 2 + (1 - 9 ) y^j 

^ 0K+1 (yi^ ((’^ + 1)”) > 2 / 2 ) + (1 - (y[Y (^(n+ 1)") , 2 / 2 ) • 

One can use the same technique to show that monotonicity is preserved. An identical argument can be carried 
out for the GMWB. 

Convexity and monotonicity preservation are established for a stochastic volatility model in [4]. For 
the case of general parabolic equations, convexity preservation is established in [16]. This result is further 
generalized to parabolic integro-differential equations, arising from problems involving assets whose returns 
follow jump-diffusion processes [5]. 

4.5 Existence of an optimal bang-bang control 

Once we have established that convexity and monotonicity are preserved across and between exercise times 
(i.e. §4.3 and §4.4, respectively), we need only apply our argument inductively to establish the existence of an 
optimal bang-bang control. Instead of providing a proof for the general case, we simply focus on the GLWB 
contract here. For the case of a general contract, assuming the dynamics followed by the assets preserve the 
convexity and monotonicity of the cost of funding the contract between exercise times (e.g. GBM, as in 
Appendix A), the reader can apply the same techniques to establish the existence of a bang-bang control. 

Example 4.22. Consider the GLWB. Suppose that for some n s.t. 0 ^ n < N, Example 4.21, 

is also CM. Under sufficient regularity (see Appendix A), for fixed x, Vx,n Is bounded above (satisfying 
{A2)). Since (Al) and (A2) are satisfied, we can use Example 4.14 to conclude that the supremum ofvx,n, 
for each x G 11, occurs on {0,1,2}. By Example 4.20, V~ is convex and monotone. 

By (2.10) and (2.11), V (x, A^) = 0. Since V (x, A^) is trivially CM as a function of:x., we can apply the 
above argument inductively to establish the existence of an optimal bang-bang control. 

5 Numerical Examples 

To demonstrate the bang-bang principle in practice, we implement a numerical method to solve the GLWB 
and GMWB problems and examine loss-maximizing withdrawal strategies. 

5.1 Contract pricing algorithm 

Algorithm 1 highlights the usual dynamic programming approach to pricing contracts with finitely many 
exercise times. Note that line 2 is purposely non-specific; the algorithm does not presume anything about the 
underlying dynamics of the stochastic process(es) that F is a function of, and as such does not make mention 
of a particular numerical method used to solve Vfi' given l^j. Establishing that the control is bang-bang 
for a particular contract allows us to replace A„ appearing on line 4 with a finite subset of itself. 

5.2 Numerical method 

The numerical method discussed here applies to both GLWB and GMWB contracts. Each contract is 
originally posed on il = [0,oo)^. We employ Algorithm 1 but instead approximate the solution using a finite 
difference method on the truncated domain [0,a:™“] x [0,a:™’^’^]. As such, since fx.n (A) will not necessarily 
land on a mesh node, linear interpolation is used to approximate Ififi (fx,n (A)) on line 4. A local optimization 
problem is solved for each point on the finite difference grid. Details of the numerical scheme can be found in 
[ 1 , 10 ]. 
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Data: payoff at the expiry, Vn = 

Result: price of the contract at time zero, Vb = 

1 for n ^ N — 1 to 0 do 

2 use to determine 

3 for X € n do 

4 I V- (x) EE SUP;^gA„ Vn (fx.n (A)) + /x,n (A) 

5 end 

6 end 

Algorithm 1: Dynamic programming for pricing contracts with finitely many exercise times. 

Between exercise times, the cost of funding each contract satisfies one of (2.8) or (2.16). Corresponding 
to line 2 of Algorithm 1, we determine from using an implicit finite difference discretization. No 
additional boundary condition is needed at Xi = 0 or 0:2 = 0 ((2.8) and (2.16) hold along dfl x +i))- 

The same is true of X 2 = » 0. At xi = » 0, we impose 

Vixr^X2,t)=g{t)xr^ (5.1) 

for some function g differentiable everywhere but possibly at the exercise times n. Substituting the above 
into (2.8) or (2.16) yields an ordinary differential equation which is solved numerically alongside the rest of 
the domain. Errors introduced by the above approximations are small in the region of interest, as verified by 
numerical experiments. 

Remark 5.1. Since we advance the numerical solution from n~ to (n — 1)’*’ using a convergent method, 
the numerical solution converges pointwise to a solution V that is convexity and monotonicity preserving. 
Although it is possible to show—for special cases—that convexity and monotonicity are preserved for finite 
mesh sizes, this is not necessarily true unconditionally. 

Remark 5.2. Although we have shown that an optimal bang-bang control exists for the GLWB problem, we 
do not replace A„ with {0,1,2} on line 4 of Algorithm 1 when computing the cost to fund a GLWB in §5.3.1 
so as to demonstrate that our numerical method, having preserved convexity and monotonicity, selects an 
optimal bang-bang control. For both GLWB and GMWB, We assume that nothing is known about Ux,n o,nd 
hence form a partition 

Ai < A2 < • ■ ■ < Ap 

of the admissible set and perform a linear search 

5.3 Results 

5.3.1 Guaranteed Lifelong Withdrawal Benefits. Figure 5.1 shows withdrawal strategies for the holder 
under the parameters in Table 5.1 on the first four contract anniversaries. We can clearly see that the 
optimal control is bang-bang from the figures. At any point (x,n), we see that the holder performs one 
of nonwithdrawal, withdrawal at exactly the contract rate, or a full surrender (despite being afforded the 
opportunity to withdraw any amount between nonwithdrawal and a full surrender). 

When the withdrawal benefit is much larger than the investment account, the optimal strategy is 
withdrawal at the contract rate (the guarantee is in the money). Conversely, when the investment account 
is much larger than the withdrawal benefit, the optimal strategy is surrender (the guarantee is out of the 
money), save for when the holder is anticipating the triennial ratchet (times n = 2 and n = 3). Otherwise, 
the optimal strategy includes nonwithdrawal (to receive a bonus) or withdrawal at the contract rate. Note 
that the strategy is constant along any straight line through the origin since the solution is homogeneous of 
order one in x, as discussed in [10]. 
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Withdrawal account ( 22 ) Withdrawal account ( 22 ) 


Figure 5.1: Optimal control for the GLWB for data in Table 5.1. As predicted, there exists an optimal 
control consisting only of nonwithdrawal, withdrawal at the contract rate, and a full surrender. 


□ No withdrawal □ Withdrawal at the contract rate Full surrender 
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Table 5.1: GLWB parameters. 


Parameter 


Value 

Volatility 

a 

0.20 

Risk-free rate 

r 

0.04 

Hedging fee 

a 

0.015 

Contract rate 

s 

0.05 

Bonus rate 

p 

0.06 

Expiry 

N 

57 

Initial investment 

Wo 

100 

Initial age at time zero 


65 

Mortality data 


[20] 

Ratchets 


Triennial 

Withdrawals 


Annual 


Anniversary n 

Penalty k„ 

1 

0.03 

2 

0.02 

3 

0.01 

^4 

0.00 


5.3.2 Guaranteed minimum withdrawal benefit. For the GMWB, (Cl) is violated. In particular, for 
Kn > 0, the function (A) is concave as a function of x. However, when = 0 or G = 0 (G = 0 is 
considered in [15]), the function (A) (see (2.7)) is linear in x, and hence the convexity of V~ can be 
guaranteed given CM. In this case, it is possible to use the same machinery as was used in the CLWB 
case to arrive at a bang-bang principle (see Theorem 4.12). The case of = 0 corresponds to zero surrender 
charges at the nth anniversary, while G = 0 corresponds to enforcing that all withdrawals (regardless of size) 
be charged at the penalty rate. 

Now, consider the data in Table 5.2. Since = 0 for all n ^ 7, the convexity of H in x is preserved 
for all t G (6,iV]. However, since Kg > 0, the convexity is violated as t —)■ 6“. Figure 5.2 demonstrates this 
preservation and violation of convexity. As a consequence, V will not necessarily be convex in x as t — >■ 5"*", 
and the contract fails to satisfy the bang-bang principle at each anniversary date n < 5. 

Note that for X 2 > 0, the conditions Xx 2 S [0,GAa:2] and Xx 2 G {G A X 2 ,X 2 ] appearing in (2.7) are 
equivalent to A G [0, G/a:2 A 1] and A G {Gjx 2 A 1,1], respectively. Assuming that is CM and taking 

(x) = {Pi, P2} with Pi = [0, G/x 2 a 1] and P2 = [G/a;2 A 1,1] yields that there exists an optimal control 
taking on one of the values in (0, G/a;2,1} at any point (x,n) with X 2 > 0. These three actions correspond 
to nonwithdrawal, withdrawing the predetermined amount G, or performing a full surrender. This is verified 
by Figure 5.3, which shows withdrawal strategies under the parameters in Table 5.2 at times n = 6 and 
n = 7. As predicted, along any line X 2 = const., the optimal control takes on one of a finite number of values. 
Since at n = 6, k„ > 0, we see that the holder is more hesitant to surrender the contract whenever xi ^ X 2 
(compare with the same region at n = 7). Control figures for CMWBs not satisfying the bang-bang principle 
can be seen in the numerical results in [9, 7]. 

Remark 5.3. Consider a GMWB with k„ = 0 for all withdrawal times n. As suggested by the above, this 
contract satisfies the bang-bang principle (in particular, Theorem 4.12 is satisfied) everywhere. However, the 
GMWB does not necessarily yield an optimal bang-bang control since V(f (x.) depends on X 2 (Corollary 4.13 
is not satisfied). For example, consider an optimal control for the GMWB taking on the value Gfx 2 at each 
X with X 2 > 0. Such a control’s range is a superset of (0,oo) (not a finite set). However, in this case, the 
bang-bang principle guarantees that for fixed X 2 , a finite subset of the admissible set need only be considered in 
the corresponding optimization problem. Computationally, this is just as desirable as the case of an optimal 
bang-bang control. 
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Figure 5.2: V (x, i) for fixed x\ = 100 under the data in Table 5.2. Points where V (x,n ) = F 
correspond to nonwithdrawal. To the left of these points, the holder performs withdrawal (see Figure 5.3^. 



to t ^ Q . 



t^l-. 


Table 5.2: GMWB parameters [8]. 


Parameter 


Value 

Anniversary n 

Penalty 

Volatility 

a 

0.15 

1 

0.08 

Risk-free rate 

r 

0.05 

2 

0.07 

Hedging fee 

a 

0.01 

3 

0.06 

Contract rate 

G 

10 

4 

0.05 

Expiry 

N 

10 

5 

0.04 

Initial investment 

Wo 

100 

6 

0.03 

Withdrawals 


Annual 

^ 7 

0.00 


Figure 5.3: Optimal control Ax scaled by X2 for the data in Table 5.2. 


□ No withdrawal □ Withdrawal of min (x2, G) ^ Full surrender 
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6 Conclusion 


Although it is commonplace in the insurance literature to assume the existence of optimal bang-bang controls, 
there does not appear to be a rigorous statement of this result. We have rigorously derived sufficient conditions 
which guarantee the existence of optimal bang-bang controls for GMxB guarantees. 

These conditions require that the contract features be such that the solution to the optimal control can be 
formulated as maximizing a convex objective function, and that the underlying stochastic process assumed 
for the risky assets preserves convexity and monotonicity. 

These conditions are non-trivial, in that the conditions are satisfied for the GLWB contract but not for the 
GMWB contract with typical contract parameters. From a practical point of view, the existence of optimal 
bang-bang controls allows for the use of very efficient numerical methods. 

Although we have focused specifically on the application of our results to GMxB guarantees, the reader 
will have no difficulty in applying the sufficient conditions to other optimal control problems in finance. We 
believe that we can also use an approach similar to that used here to establish the existence of optimal 
bang-bang controls for general impulse control problems. In the impulse control case, these conditions 
will require that the intervention operator have a particular form and that the stochastic process (without 
intervention) preserve convexity and monotonicity. We leave this generalization for future work. 


A Preservation of convexity and monotonicity 

In this appendix, we establish the convexity and monotonicity of a contract whose payoff is CM and written on 
assets whose returns follow (multidimensional, possibly correlated) GBM. We do so by considering the PDF 
satisfied by V and the fundamental solution corresponding to the operator appearing in the log-transformed 
version of this PDF. Considering the log-transformed PDF allows us to eliminate the parabolic degeneracy at 
the boundaries and to argue that the resulting fundamental solution for the log-transformed operator should 
be of the form P (y, y', r, r') = P (y - y', r, r')- 

We begin by describing some of the notation used in this appendix: 

• Let n = rii X 172 where = (0, oo)"* and 172 is a convex subset of a partially ordered vector space A 
over K with order 17 can thus be considered as a convex subset of the vector space A = M"* x A 
over R. 

• We write an element of ft in the form (x, x^+i) = {xi, ..., Xm+i) with x € Di and x^+i & in 
order to distinguish between elements of f7i and fl 2 - 

• The partial order we consider on A is simply inherited from the orders ^ on R™ (Remark 4.8) and 

Specifically, (x,a;„+i) if and only if x x' and x^+i a^m+i- 

Suppose V satisfies 



dtV + CV -1- w = 0 on 17 X (7„, 7„+i) 

(A.l) 

and 

where 

V (x,Xm+l,tn+l) = On 17 

^ m m 

(A.2) 


(A.3) 


^,j=l i=i 

In the above, uj = lo (x,t). We will, for the remainder of this appendix, assume the following: 


(Dl) aij = Gij (t), bi = bi (t), and c = c{t) (i.e. the functions bi and c are independent of (x,a;m_|_i)). 
(D2) Y^T,j=i^i, 3 ^xidxj is uniformly elliptic. 
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Example A.l. For the GLWB guarantee, C is given in (2.9) and uj = M. {t) xi. 

Remark A.2. We say V satisfies (A.l) if V is twice differentiable in (the components of) x and once 
differentiable in t on Q.X (t„,l„+i),'* continuous on x {tn,tn+i],^ and satisfies (A.l) pointwise. 

We now describe the log-transformed problem. For ease of notation, let 

,..., e^™) a'^ (r) = (t„+i - r) 

F {y,ym+i) = F{e^,ym+i) bl (r) = 6* (t„+i - r) 

{y,ym+l,T)=Uj{ey,ym+l,tn+l-T) c' (t) = C (t„+i - r) 

and n' = K™ x n 2 - Let E be a solution of the Cauchy problem (A.l) and (A.2). Let 

u (y, 2/m+i, r) = y (e^, ym+i,tn+i - r) 
and A = tn+i — tn- Then, u satisfies 

C'u — drU -I- w' = 0 on n' X (0, A) (A.4) 

and 

u{y,ym+i,o) = {y,ym+i) (A.s) 

where 

- m m 

i,j=i i=i 

Note that (D2) implies that £' is uniformly elliptic. 

In order to guarantee that a solution u to the log-transformed Cauchy problem (A.4) and (A.5) exists, 
and is unique, sufficient regularity must be imposed on ip', and w'. We summarize below. 

(El) For each ym+i, y f' (y)2/m-i-i) is continuous on K™. 

(E2) The coefficients of C are sufficiently regular. 

(E3) For each j/^+i e ^ 2 , (y,T) w' {y,ym+i,T) is sufficiently regular. 

(E4) u satisfies a growth condition as |a;| —>■ oo. 

For an accurate detailing of the required regularity, see [11, Chap. 1: Thms. 12 and 16]. 

When (D2) and (El) — (E4) are satisfied, the solution u can be written as 


u{y,ym+i,T) 


L (y, y', r, 0) ip' (y', ym+i) dy' 

+ f [ r(y,y',r,r')a;'(y',2/,„+i,T')dy'dr' on X (0,A) (A.6) 

Jo Jm’^ 


where F is the fundamental solution for C (whose construction was first detailed by [17]). We first note that 
(El) follows immediately if tp is convex, as shown below. 

Lemma A.S. If ip is convex w.r.t. the order (see Definition 4.3^, then for all x^+i, y ^ (y; Xm+i) 

is continuous on M™. 

^i.e. S (ft X 

®i.e. y\Qx{t„,tn+l] S C (C X {tn,tn+l]). 
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Proof. We have assumed that tp = tp {k, Xm+i) is convex w.r.t. on From this it follows that for all 
Xm+i S ^2) is convex in x on fli = [0, oo)™ w.r.t. to the order ^ on K"*. This in turn yields that for all 
Xm+i € ^2, is continuous in x on fli. Therefore, ip' = (p' (y; Xm+i) is continuous in y on K”^. □ 

Theorem A. 4. Suppose (Dl), (D2) and (E2) — (E4). Suppose that ip is CM w.r.t. the order (see 
Definition 4.3 and Lemma A.7). Suppose further that for all t € ^ is CM in (x,a;m+i) on fl w.r.t. 

the order Then, for all t S (tn,tn+i], V is CM in (x.,Xm+i) on LI w.r.t. the order In particular, 
V+ is CM. 

Proof, r appearing in (A.6) depends on y' and y through y' — y alone since by (Dl), a'^, 6' and c' are 
independent of the spatial variables [11, Chap. 9: Thm. 1]. Therefore 


u(y,2/m+i,T) = / r(y'-y,r,0)</5'(y',j/„+i)dy' 

+ f f r(y'-y,T,r')w'(y',ym+i,r')dy'd'r'onK"* X (0,A). 

Jo Jm.^ 

Let logx = (logxi,... ,logXm)- From the above, whenever Xi > 0 for all i ^ m, 

V {x,Xm+i,t) = / T(y'-logTi,tn+i-t,0)p(ey ,Xm+i)dy' 

+ / r(y'- logx,t„+i - t,r')a; (e5"',x„+i,t„+i - t') dy'dr' on D X . 

Jo JR^ ^ ' 

Denote by x o x' = {x\x’i ,..., Xmx’m) the elementwise product of x and x'. The substitution y' = log (x o x') 
into the above yields 


POO POO ^ 

V{x,Xm+i,t)= / ■•■/ r(logx',t„+i-t,0)(^(xox',a;m+i) ==— ydx' 

Jo Jo lli^i 

pOO pOO ^ 

+ / / •■•/ r(logx',t„+i -t,T')w(xox',a:m+i,tn+i- t') — ^dx'dr' on D x (t„, t„+i) . 

Jo Jo Jo lli^i 


Since T is > 0 [11, Chap. 2: Thm. 11] (a related, arguably more general result is given in [12, Chap. IV: 
Prop, l-llj), from the convexity and monotonicity of and ui, it follows immediately that V„ (p^,Xm +i)t) 
is CM on D for any t G (t„,t„+i). 

Remark A. 5. IFe can extend our construction to Di x LI 2 by taking limits of V up to the boundary. Since 
the codomain of V is Hausdorff, this extension is unique. 

□ 


B Commutativity of union and supremum 

Let r be a poset with order ^ satisfying the least-upper-bound property. All supremums are taken w.r.t. T. 
Lemma B.l. Let S = be an indexed family of nonempty subsets ofT. Let S = Uae^ (ind 

U = {supScl^g^. 

Then, sup S = sup U whenever S is bounded above. 
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Proof. Suppose A is empty. Then both S and U are empty, and hence the expressions agree. 

Suppose A is nonempty and that S is bounded above. Since S is bounded above, its supremum u must 
occur in T. For each a, m is an upper bound of Sa, and since Sa is a nonempty subset of T, sup Sa = Ua 
for some € T. Thus, U = C T. Since Ua ^ u for each a, u is an upper bound of U. Since A is 

nonempty, U is nonempty and hence U has a least upper bound u' G T with u' ^ u. Let x G S. Then x G Sp 
for some l3, and hence x Ujs u' so that u' is an upper bound of S. Since sup S = u, u u' and hence 
u' = u. □ 
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